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ABSTRACT 

In this paper we discuss the ionization structure of hydrogen and 
helium in dusty nebulae. The general equations of Paper I (Petrosian 
and Dana 1075 ) fo r pure hydrogen nebulae have been modified to include 
helium. We first present an approximate analytic solution for Stromgren 
radii of hydrogen and helium in the absence of dust. We then show that 
these results, with simple modifications, are also applicable n dusty 
nebulae where the effective absorption cross section of dust grains 
varies slowly with frequency in the 1000 to 20oA range. No analytic 
solutions are possible if this cross section varies rapidly with fre- 
quency. In this case, however, we have derived simple coupled differen- 
tial equations which can easily be solved numerically. We present 
approximate analytic expressions for evaluation of the variation of the 
fraction of ionizing radiation absorbed by dust and the ratio of the 
volume emission measures of He II to H II regions with the spectrum of 
the ionizing source, helium abundance and absorption properties of dust. 
The effects of dust on the He III zone are discussed in the Appendix, 

As in Paper I our results are restricted to spherically symmetric nebulae, 
but non-uniform gas and dust distributions and clumpiness can be taken 
into account by our general results. 

_ 
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I. INTRODUCTION 


In Paper I we discussed various approximate analytic solutions of 
radiative transfer equations in nebulae contained hydrogen and dust. 

In this paper we consider dusty nebulae containing hydrogen and helium. 

Our aim is to obtain analytic solution or simple differential equations 
which could be solved numerically easily. 

In section II we present modification of the general equations in 
Paper I so that they would be applicable to the present case. As shown 
here, with very few simple approximations (approximation usually made even in 
detailed numerical calculations such as the on-the-spot approximation' 
the problem of transfer of ionizing photons can be reduced to two simple 
coupled differential equations. 'Ill esc are equations 10) and ' 17 ). 

These equations form the basis for the rest of this paper. In 
section III we reconsider the pure hydrogen nebulae with and without 
dust. This allows us to simplify the treatment of the II plus He case 
which is discussed in section IV. Here we show that in absence of dust 
one can describe the ionization of H and He analytically. To our knowledge 
no such solution has been presented before. The generalization of these 
results in dusty nebulae is simple and is presented in section IVb. In 
the concluding section (V) we apply these results to two observationally 
interesting parameters: the fraction of ionizing photons absorbed br 

gas (or dust'' and the ratio of He to H volume emission measures. We 
give simple foimulae for calculation of these quantities for the general 
problem of dusty nebulae. 


II. GENERAL EQUATIONS QF TRANSFER OF IONIZING PHOTONS 


'Hie general equations of rndintivc transfer and Ionization 
equilibrium described in section II. 1 of Paper I can be applied 
to nebulae with hydrogen, helium and dust with following 
modifications . 

In equations (I. I 1 ) and (1.3) (all figures and equations of Paper I 
will be identified with the prefix i) we must include absorption and 
emission coefficients of helium, 


K v , no “ i l -y)vn® no (v) 


Ur j 


V; 


He 


“ ynn eK,v + Cl * 


cr 




1 > 


where n is the total hydrogen (H I and H II ) density, 


n 


n( He II) 
y “ n(lle) ~ 


n = ( x + y Y ) n 
e 


( ;>> 


^ is fraction of photons from recombination to excited states of helium 

capable of ionizing hydrogen ( varies between 0.8 and O. 96 , cf. e.g. Mathis 

/ ( 2 ) / 

I 97 I), and and or 1 are the recombination coefficients to the 

(i)' 

ground and excited states of helium («' t= a' + or"' ). Similarly 
equation ( 1 . 6 ) must be supplemented by the ionization equilibrium of 
helium: 


4* 
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J K „dv = xnn (X , 4jt 
v Vr e 
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v 0 


1.8v , 


J K „ dy = yYnn a . : 3' 
v vjHe J e 


'•’0 0 

It should also be noted that in the presence of helium the upper 
limit of integrals over freouency will be 4 m instead of co , because 


photons with v > 4^ will be absorbed in the inner zone where He will 
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bo doubly ionized. For H II rogions this zone is in general negligible. 
For nebulae with central sources omitting considerable numbers of Ho II 
ionizing photons our results will apply to the rogions beyond this zone 
with minor modifications (cf. Hummer and Seaton 1964 ), The effect of 
dust on the He III zone is discussed in the Appendix, 

As in Paper I , we shall be dealing with intensities and fluxes 
integrated over frequency. We must however distinguish between photons 
capable of ionizing only hydrogen and those capable of ionizing both 
hydrogen and helium. We therefore define net fluxes crossing spherical 
shells 



and similar expressions for intensities I and J. 

If Sp photons were absorbed by helium alone, the equations for 
H and He ionizing photons would be decoupled and the results of Paper I 
would apply to S-^ and Sp photons separately. However, because H 
and He compete for Sp photons, the solution of the problem is more 
complicated. Instead of equation (1.9) we must define separate average 
hydrogen cross sections for v < 1.8v Q and v > l,8v Q photons: 


r‘ 8v ° 3 r 4v ° , 

1 °H,1 “ I ctq(v/vq) J(v)dv > J 2 CT H,2 = J Oq(v/Vq) J(v)dv 
v 0 1,8v 0 {5) 


so that the total average hydrogen cross section is 


4 




* D 

As in Paper I, we can define separate cross sections a and a for 

X X 

•X* 

the stellar and diffuse radiation. For cr x the average intensities 

*X > 'X* 

J 0 and can be replaced by and Sp respectively. 

We shall use the on-the-spot approximation and set the dust albedo 
<d(v) - 0* As we have shown in Paper I, the effects of absorption by dust 
of diffuse radiation and scattering by dust can be accounted for by 
defining effective dust absorption optical depths. The relationship 
between the effective and the true absorption optical depths can b - 
found in figures 1 .^ and 1 , 7 , 

With these approximations the diffuse radiation satisfies the 
relation 

K He^ 1+ P ^ ItJ 2 = Yynn e a l ' K H,1 4nJ l = m e a i + Cx Y y nn e a ^ > (9) 
where 

p = k h ; 2 /K He “ x )o- h (!.8vo)/[ y ( 3- - ( ls ^) 
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is the fraction of recombination photons to ground state of helium 
absorbed by hydrogen. Substitution of equation ( 9 ) ln equation (3' 
for ionization equilibrium gives 

(1-x) - A^(x + yY)[x - yY((J^+ £->)] > (1-y) = A^y(x+ yY) (1+ £,,) (11) 

whero we have defined 

V = Ya^ /a^ , C 2 85 a l P/[ a ^ ( 1+ 0 )] > d' 1 ) 


and 


A, . W 2 n a (2) /(s*4) 


A n = 


* (o)* * ( 'A 

Vrf*'"’ ''<YV‘ > 


. (13) 


Substitution of these in the equations governing transfer of stellar photons 
(dS*/dr . -K tot;1 S*) gives 


dS-, 


3 l /dr = - K d,l S l - ^r 2 nn e a (2) [x-Yy(C 1 + C 2 ) ]4,l4 / 4 S * ' 

(lb) 

dS*/dr = -Kj 2 S* - Wnn e a^^(l+ C 2 )yY+ [ x " yy Ui H ' £ 2 ^ ]<4, 2 S 2 /(7 H S } 


Addition of these two equations gives 

dS*/dr = -K d S* “ ^rtr‘*nn e 0 f^[x+ yY(l - £^)] , (19) 

y ' 

In most treatments of the problem this fraction (and consequently the 
quantity in eq. [12]) is set equal to zero. If -his was the case 

the absorption of S 0 photons by hydrogen could also be neglected. As 
we shall see below, neglecting p or with respect to unity will 

cause up to 30 percent under-estimation of Yq or a similar over- 
estimation of the relative He abundance Y. 
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which for = 1 is identical to equation ( I 
has a contribution from helium (note also that 
on y) . 


,2‘j) except that horp n p 

■H* # 

tC d A d ^ s 1 + <* Y doponds 


From here on we shall bo dealing only with stellar photons, there- 
fore we shall eliminate the * notation in what follows. With this 
modification and with elimination of S, and S 0 in favor of Y c S,,/S 
we find the general expressions for stellar photons are 


dS/dT = -(1+ cey)S - f(r) , 


dy/dT = -otY(l-Y) + (1- y)(PY-Y , )(1-py)“ 1 f(r)/s 


( 16 ) 


where 

f(r) = [x+ yY(l - C 1 )]4nr 2 nn o a (2) /K d ^ 1 , d T = K^dr 

yY(i+ c ? ) 

y' _ _£•_ 

x+ yY(l - Qj) 

Another approximation from Paper I which we employ here also is to 
ignore the (1 - x) and (1-y) terms in equations (1^) to (I 7 ) with respect 
to unity. 11113 is true everywhere (in particular for large dust optical 
depths) except very near the ionization fronts (cf. fig. 1.2). This 
amounts to setting x = y = 1 in equations (14) to (I 7 ). 


} 

( 17 ) 
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III. PURE HYDROGEN NEBULA 


(a) Nebulae Without Du at 

In the absenco of dust and helium (K^ = y = Y = 0) equations (11) 
and (l^), with the approximations stated in section I , give 

a % - 1 - 5 * f 1 • - (* + 35) - [( l + 5^)“ - x ] 1/2 < lf ’> 

where is defined in equation ( 13 ) t 

r 

?' “ - 4nr 2 n 2 a ( 2 ) /S 0 , g(r) = J ?'dr (I 9 ) 

r 0 

and 


S q - | 4nr 2 n 2 a^ ^dr (20) 

■0 

is the total number of ionizing photons emitted by the central source(s). 
Here r^ and r^ are the inner and the outer radii of the ionized 
region respectively; ^(r^) = 1 and in general 1*0 « These are 

identical to equations ( 1 , 21 ) and ( 1 , 22 ). 

(b) Nebulae With Dust 

In this case K / 0. It can be shown that results identical to 
a 

equations (18) to (20) are obtained if we replace S and §' by 

S = Se T and 'f' = f»'e T (21) 

T 

and if we include e in the integrand of equation (20); dT = n j 0 -,dr. 

d d 

The results thus obtained are identical to that of equations (l.2u) to 

(1.32). 
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IV. DUS TV NEBULAE CONTAINING HYDROGEN AND HELIUM 


Since to our knowledge there are no approximate analytic solutions 
for nebulae with H and Ho (but without dust) wo consider this case 
first. 

(a) Nebulae Without Dust 

With => 0, equations (16) and (17) reduce to 

S/S 0 - (1-5) , dY/dS = - I (1 ? > (“) 


where and § are now 


r l 

S Q -f [1+Y][1+ Y(1 - Ci)3 ^r 2 n 2 a (2) dr , 

*0 

Z = J [1+ Y][l+ Y(1 - g x )] lter 2 n 2 a^ 2 ^dr/S 0 . 


(23) 


Figure 1 shows the position of r^, r g and r^, the Strb'mgren radii of 
++ + “♦* 

He , He and H , respectively. Within our approximation r ^ g r^. In 
cases where r n < r^ the parameters Y and Y must be set equal to 
zero for r^ > r > r^. 

Hie parameters (3 and Y 7 in equation ( 22 ) vary throughout the 
nebula. Variation of (3 is due to change in the spectrum of ionizing radia- 
tion, However, since the absorption cross section of hydrogen decreases 

-3 3 

rapidly with frequency (j «v ), p 1 - (1/1.8) ~ 0.8 and changes 

by a few percent for a variety of plausible spectra (cf. table 1). We 
therefore neglect variation of ( 3 . Ihe parameter Y', on the other hand. 
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varies because of tho variation of Qg which is primarily duo to the 
change of the ratio (l-x)/(l-y) in equation (lo). In general, if y < Y, 
Qr, « 1 and it can bo ignored. However, for Yq > Y, Qg > 0.1 in tho 
inner regions, and near tho edge of tho nebula whoro y - 1, -* 

u 

~ 0.6 (for electron temperature of 10^ °K) , As wo Shull see 
bolow, oven in this case tho variation of Qg iB negligible. Thus wo 
shall assume that Y' is a constant. 

Ionizing Flux 

With tho above assumptions, equation (22) is readily solvod: 



„ /l-Y t 1 ' 1 ' / BY 0~ Y ' \ 1 - Y ' 
\1-Yq/ 'py-y' ' 


(24) 


Figure 2 shows the variation of y with S for various values of 
Y 0 (solid lines). For PY Q = y' , Y = Y Q = constant. For fiY Q < 
Y (and therefore S£) becomes zero at 

s = S cr = s 0 (l - PYqW 1 - - Y ' ) (1 . Y 0 ) (P - 1)/(p - Y ' > 


Y' 


( 25 ) 


or at r = Yg, whore r 2 is obtained from equation ( 23 ) with ? 8 
si- S cr /®Q’ ’Hius, in these cases helium ionization stops before the 
edge of the nebula. For PYq > Y / , Y increases toward the outer edge 
and approaches unity at the edge of the ionized region (at r c r^), 
where S = 0. As is ovident, the shape of these curves is determined 
primarily by the value of Y 0 (actually by the value of PYq/y' ) . 
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Ionization Structure 


Onco tho variations of S and Y aro known tho ionization structure 
cun be calculated according to equations (11) to (13’'. Sinco 1.0 < 

r > 

+ £ rj ) £ 1,4, equation ( 11 ) can bo approximated ns ( 1 -x) A^x' 
so that ( 1 -x) is given by equation (lB) with a ^( 1 -(V) 

(cf. oq. [6 ]), A few valuoB of eT|j -j/qq “bd Cj^/cTq aro given in 
table 1. The slow variations of those quantities duo to changos of 
the spectrum throughout the nebula are neglected in the above 
treatment. 

Equation (11) can now be solved for (1-y) using the above values 
of x and assuming £ 0 = constant. This assumption is clearly Justified 
since A r , varies much more rapidly than any expected variation of £ 0 , 

Hie results of this calculation are showi on figure 3 for various values 
of Y » 0.1 and uniform nebula. Here we also show variation of , 

As evident; the assumption of constancy of is a good approximation 
for small values of Y n . Fo • 1 values of Y r ,, Co « 1 90 that 
y' m y" m l.O^Y (for electro t n . «<res of about 10^ °K, cf. Burgess 
and Seaton I 960 ). For large calu ^ varies slowly but it is 

no longer negligible compared to in. j. Therefore, neglecting as 
is commonly done, will imply Incorrect determination of the value of 
Y^/y' from observations. For example, for a given value of v and ob- 
served value of r 2^ r l f the re£ l uirec * value of Yq underestimated when 

is neglicted. In general for small r^/v^y Co negligible and it 

rarely exceeds 0.4 (cf. figs. 5 and 7 ). 

(b) Nebulae With Dust 

In this case simple analytic solutions vo the coupled differential 
equations (l 6 ) and (I 7 ) are possible only if G is zero or negligible, 


11 



i.o. only if a d (v) varies slowly wo that c d mi ct ( , ^ mi c d We 
discuss this cnso firBt, 

(i) « 0 

In this case, as in soction II. b, wo can dofino £5 and as 

in equation (21) and obtain rosults identical to those presented by 

equations (22) to (2|;) with replacement of S and Z* by s’ and tj' 

and with inclusion of o T in tho intogrand for S Q , Thus, the solid 

linos in figuro 2 give the variation of y with S. However, variation 

of y with S will be different (dashed lines). Using these results 

we have calculated variations of (1-x), (1 - y) and for uniform 

nebulae (n d « n o constant) and for a dust to gas ratio and gas column 

-4 

density such that t-^ - n j^ d r i % 1.0 and n d c d / n a 0 s e = 2,(j x 10 . 
These results nr> presented in figuro 4. Comparison of figures 3 and 4 
show that fr-.v a given size r^ of tho nebula, introduction of dust 
increases the fractional ionization of both hydrogen and helium. This 
is primarily duo to increased values of S throughout the nebula. 

(ii) a J 0 

In this case simple analytic solutions are not possible. However, 
equations (16) can readily be solved numerically. We have solved those 
equations for few values of (X, y^ and Y, As above, these results can 
be used to obta-'n the ionization structure of the nebulae. Some of these 
results are presented in figures 2 and 5 where we show the variation of 
^2 with a for various values of y Q and Y. 

In general we find the shapes of (1-x) and (1-y) curves are 
nearly independent of the details of the problem and are determined 
primarily by the vr js of r 2 and r, . Consequently, we do not present 


* 


12 


tho variation t f (1 - x) and (1 - y) fo* values of 01/0, 
we concentrate on the dependence of observable quantities on 
other parameters In the next section. 
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k 


Instead 
a and 


V. SUMMARY AND RESULTS 


We have derived approximate analytic solutions to the equations of 
radiative transfer .\nd ionization structure of nebulae containing hydro- 
gen and helium. These solutions can easily be extended to nebulae con- 
taining dust if the effective dust absorption cross section is a slow 
varying function of frequency in the 1000 to 2(X)A wavelength range. For 
dust with a widely varying cross section in this wavolongth range we have 
derived simple differential equations vhich can readily be solved numerically. 
The main approximation of our treatment is neglect of albedo of the dust 
grains. The other approximations are of minor consequence. 

Since we treat fluxes integrated over wide ranges of fiequency, we 
lose most of the information on the variation o t the spectrum of ionizing 
radiation throughout the nebulae. Consequently, we cannot calculate the 
ionization structure of heavier elements. The main result of these solu- 
tions is the ionization structure of He and H. We summarize these results 
below. 

(a) Fraction o f Ionizing Radiation Absorbed by Gas 

A useful parameter for comparison with observations is the fraction 

of stellar ionizing photons absorbed by gas (or dust). The fraction 

f . absorbed by gas is 
net 

r 1 

f net = 4 * a(2) /* (1+ Y)[l+ Y(l- ; 1 )]r 2 n 3 dr/S 0 . (26) 

r 0 

For a = 0 and ^ m 1 equation (26) reduces to 


Ik 


not 


f(T x ) + Yf(Tg)g(T 2 )/g( Tl ) 

1 + Yg(T 2 )/g(T 1 ) 


( 27 ) 


whore (cf. eqs. [ 1 . 27 ] and [1.32]) 


K ^ T 

g(T) f n 2 r # 2 exp(T / )dr / /K d , 

n o Jo 

f(T ) . ^ nW/K . 

&(T)n 0 Jo 


( 28 ) 


In these equations n Q is the value of n at r s r Q , and and 

T g are the dust optical depths for photons up to the H and Ho 

Stromgren radii r^ and r 2 , respectively: 



(29) 


The fraction f , calculated from equation ( 27 ) is nearly equal to 
net 

f ( ) ) the fraction obtained in Paper I for Y = 0. Ibis equality is 
exact for Tg/T^ = 0 and 1.0. In the Intermediate range the 
difference between the two cases is less than three percent (ior 
£ 5) because Yg(Tg)/g(T 1 ) <Y mt 0.1. 

For a 4 0 equation ( 27 ) is not valid because Sq is no longer a 
simple function of and T g . In this case f is given by 

(cf . also Mathis 1971) 


where 


f and 
S 1 


f net= ( 1 -V 0 ) f s 1 + \ % 

f_ are fractions of S-. and 
S 2 


(30) 

photons absorbed by 
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the gas • From equation (14) these are (for fa w 1) 



H (1+ Y)(l - Y' n^r^dr + j n 2 r 2 dr 

r 0 r _2 

/ (i + v)(i - v' )(^L) ” 2rVdr * J 1 " Vo 


'dr 


and 




, /n Y/Y' 1 2 2. 

+ ( X -P) TTpyJ n r dr 

, Y/Y'* 1 2 2 (1+ «)t . 

+ TTpyJ n r e ^ 


(3D 


Examination of the results presented in figure 2 will show that the 
quantity (l-Y)/(l-pY) is nearly a constant (approximately equal to one) 
except for large values of Yq where it goes to zero rapidly at the very 
edge of the nebula. The quantity in the square brackets of the expression 

for f is constant for Y <0.1 and for large values of Y~ it 
s 2 0 0 
varies slowly from 2 (at r =“ r^) to 6 (at r = r^) . Because these varia- 
tions are much slower than the variation of the remaining quantities in 

the integrands, the variations of Y can be neglected so that f <*, f(T,) 

S 1 1 

and f g = f[(l+o:)T 2 ] [ of . eq. (28); note that we have neglected 

Y - y' m Y£ 2 in comparison with unity for cases when r £ < r^ which is 
Justified since in these cases Y£ 2 < 0,03]. Equation (30) then be :omes 


f net' s(1 " Y 0 )f < T l ) + Vt (1+a)T 2] * (32) 
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The results obtained from this equation are within <d% of the results 
from the exact equation ( 26 ). 


( b) He to H Line Intensity Ratio 

Another important parameter for comparison with observation is the 
expected ratio of He to H line intensities. These ratios are easily 
obtained from the ratios of He to H volume omission measures 


R =s I r^nn e dr j\ 
*2 ' r 0 


r nn dr , 
e 


r T 2 2 // 1 

R' = I r nn dr /I 

K \ 

0 0 


r nn dr . 
e 


(33) 


If the He III zone is negligible r^ •* r Q and R 7 a 0. We shall assume 
this to be the case (cf,, however appendix). On figure 6 we plot R 
versus Yq/ y for nebulae without dust (dotted line) and for dusty nebulae 
with various values of the dust cross section parameter a. As explained 
previously (Petrosian 1973# 1974) for given Yq/Y the presence of dust with 
ct c j( v) constant (i.e. 0! a 0) increases the value of R with increasing 
values of the total optical depth (eq. [29]). For negative values 

of d (a d x > ^d ?) fewer He f° nizin g photons compared with H ionizing 
photons are absorbed by the dust and R is larger. The reverse occurs 
for positive value of a. It should be noted that R depends on 
Yq/Y and is insensitive to the value of Y. R is also fairly indepen- 
dent of non-uniformities or inhomogeneities in the gas and dust distribu- 
tions as long as the dust to gas ratio is constant. However R changes 

2 

for dust and gas distributions which are not the same. 


^Note that because of the approximations x=l for r^r^ and y=sl for r g r^ 
the ratio R s l. However, for large values of Y 0 (see figs. 3 and 4) where 
r 2" r l an( * (•** " y) < (l-x) the line intensity ratio R could be slightly 


larger than one. 



Thor© are no exact analytic expressions for R for th* genera. 


caso. However, the following recipe sooms to give fairly accurate results 
(cf. also Mathis I97I) • 

In the absence of dust, equations (22), (23), (25) and (33) give 
for PY 0 /y' < 1. 

R - 5 2 - h (V Q ) e l-(l-pY 0 /Y') (1 " Y/)/(P " Y/) (l-Y 0 ) (P ‘ lVrp ” Y/) ( 3*0 


and R *» 1 for PYq/Y'’ J» 1, 

In general in the presence of dust with arbitrary values of a it 
can be shown that the ratio R is given by 


"a-h^net ■ <»> 

For Of = 0 there exists an exact analytic solution since in this case 

?2 = f('Tg) '$2’ f net * f ^ T l' 1 and ^2 “ h ^ Y 0^ ( cf * thcj discuSRions in 
part (a) above and Section III-b-i), so that 


R a=0 = f(T 2 )h(Y 0 )/f(T l ) * (36) 

For a 0 there is no simple analytic expression for in 

terms of Yq* However, from equation ( 31) we can write = Y o^s 

where ^X -1 ^ is the average value of the inverse of the quantity in the 

square brackets. As mentioned before, f Q « f[(l+a)T ? ] and for small 

b 2 

values of Yj X 1, so that for PYq/Y' « 1 (i.e. for R « 1) we 

have 


“o- V[( 1+ “> T 2] /y ' f net • (37) 

Furthermore, for small ^Yq/y' equation ( 3*0 gives R =s Yq/Y / so that 
we can write R^, » YqA' if we define 
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Y i " V 0 f[(l+a)T 2 ]/f 


' 3 ^ 


This relationship, when generalized to all values of Y^, suggests 
that for dusty nebulae with arbitrary a 

B « = hCY 0> • '39' 

Comparison of R^ from equations (38) and (39) with our calculations 
presented on figure 6, shows that for low values of optical depth 
(T x ^ 2) and for (y^/Y) c* 3> these two results agree to within a 
few percent. Ten to fifty percent difference is obtained for large 
values of a and for y /Y' > 3. 

Finally examination of figures 5 and 6 shows that the quantity g Q 
depends primarily on R. Ibis is shown on figure 7 . This simplifies 
comparison with observations since which enters as a parameter in 

our discussion, can be determined purely observationally . Application 
of these results to observations will be discussed in Paper III of this 
series. 
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APPENDIX 


If the central source omits Ho II ionizing photons ( v > 4 Vj ) 


S„ «= f S(v)dv , v' e svs A , 1 ' 

\ 

then there will be an inner region of doubly ionized helium (r^ > r > r . 
in figure 1). Assuming on-the-spot approximation for the He II recom- 
bination photons we find 


2 (o\" 

dS,/dr = -k . .8- - 4jtr Yy'nn cr 1 *' 

3 d,3 3 e 


(A. 2) 


where k . is the average dust absorption coefficient for S„ photons, 

Q* 3 <5 

(r*)" 

yf _ n (ue IIl)/n(He) is the fraction of doubly ionized helium and a 
is the recombination coefficient to the excited states of He II. In 
writing equation (A. 2) we have neglected absorption coefficient of H and 


He I with respect to that of He II at these frequencies. This is 
Justified since for Y =» 0.1, ^^He II ** (l-x)/(l-y / ) « 1 

and K He I(3 A„ e „ «, 0.2 (1 -y)/(l-jr') « 1. 

Clearly the effect of dust is to absorb some of the photons and 

reduce the size of the He III zone as in the case of pure hydrogen nebulae. 
The details of this are obtained from equations similar to equations (18) 
to (21) with replacement of S by S^, r^ by r(,, x by y' and t 
by an optical depth appropriate for photons. 

There will also be absorption of and photons in this region. 

Absorption of S 2 photons by He I will be negligible. However because of 

recombination of hydrogen (whose total number in this region is equal to the 

2 (2) 

integral from r n to rl of 4irr xnn Cr ; ) there will be some absorption of 
0 <- © 
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and S,j photons by hydrogen. On the other hand, there will bo some 
emission of such photons from recombination of Ho III. TTie number of 


such photons have boon calculated by Hummer and Seaton (1964), For 

relative helium abundance Y M 0.1 this number is approximately equal 

to the number of hydrogon recombinations. Consequently, absorption by 

hydrogen can also be neglected. The remaining absorption by dust of S-^ 

and Sg photons can thon easily be taken into account once the effective 

dust absorption coefficients K . , and K . ,, are known. Thus, if S, 

ci j -L a, i. j *- * 

and S n n are total number of S n and S 0 photons emitted by the 

Cm y V X U 

central source, then at r(, their numbers will be reduced by factors of 

i rK j 


[ r r 2 j ( r r 5 ) 

exp j I K d ^drj and exp ^ I K d 0 dr| , respectively. Substitution of 

(•*0 > ) (Jr Q ) 


these reduced fluxes for S and replacement of r by r' n in the 

X j U v u 

equations of Sections 1 and 111 will give the correct result. 

For Y / ^ 0 the ratio R 7 ( oq , [ 33] ) which is related to the 
He II to H line intensity ratio will no longer be zero. If we define 


o/ - Xflii . i 

a d,l 


and 


/2 

r 2 - J *u,l dr# 


(A. 3) 


then it can be shown [eqs. (26 ), (28) and (29 )] that 


R' 


Yq f[(l+ of ) Tg] 

“ f < T l> 


l A. 4) 


1 ( °) 3 

where we have used the approximation S^ =* (1+ Y)4ircr ' g( t^) 1 and 

assumed r^ « r-^ so that n ~ constant for r < r^,. This is similar 
to equation (37) for 
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TABLE 1 


AVERAGE CROSS SECTIONS AND STELLAR PARAMETERS 


T(°K) + 


8 

C H,l /c 0 

c \ lo /a 0 

90300. 

0.500 

0.837 

0.452 

0.688 

62200. 

0.300 

0.808 

0.487 

0.832 

50000. 

0.194 

0.793 

O.5U 

0.900 

4oooo . 

0.108 

0.787 

0.555 

0.959 

40000. ++ 

0.271 

0.825 

0.519 

0.797 

37500. 

0.0887 

0.786 

0.564 

0.979 

30900. 

0.0445 

0.785 

0.602 

1.023 

30000. ++ 

0.0039 

0.664 

0.628 

1.138 

n,2* 

0.553 

0.916 

O.5O7 

0.410 

n= 4* 

0.169 

0.855 

0.585 

0,754 

* 

Averaged 

over power law distribution 

V s 0 !v 0 A,) ”- 

+ 

Averaged 

over black body distribution B (T), 

V 

^Averaged 

over model atmosphere (Auer 

and Mihalas I972) 


FIGURE CAPTIONS 


Fig. 1 
Fig. 2 


Fig. 3 
Fig. 4 
Fig. 


. Tno StrHmgron radii of Ho , He and H . rj is tho inner 
boundary of the nebula, 

. y, tho number ratio of helium ionizing photons to total 

ionizing photons, as a function of the normalized stellar 

flux s / s q* Y 0 * 0.^0, 0,194 and 0,044^ correspond to 

black body stollar temperotures of 9°>3 00., ^0,000. and 

30,900. °K respectively. The solid lines are for nobulao 

without dust (t^ « 0,0), Four different dust cross-Bections 

were used: a = -1 (dash-dot line), a ** 0 (short dash), 

Q a 1 (dot) and Ct o [) (long dash) where 1+ a 3 j . 0 /a . , . 

a, c- a, 1 

Two cases where w 2.0 are plotted and labeled. All 
other cases are for * 1.0. All curveB are for a uniform 
dust and gas distribution. 

. The neutrrl fraction of hydrogen (solid lines) and helium 
(dashed lines) and (dash-dot) for uniform dustless 
nebulae as a function of the normalized radius r/r^. r^ 
is the radius of the Stromgren sphere of hydrogen. The open 
circles denote the volume average of Y ■ 0.1, 

. The neutral fraction of hydrogen and helium and Q 0 as a 
function of the optical depth of dust T for uniform dusty 
nebulae with a constant dust cross section (a = 0) . s= 
n d a d r l * and Y = 0.1 for all cases. The symbols are the 

same as in Fi^, 3. 

1, The variation of average (over volume) with a for Y = 0.1. 

For = 0,30, we show the variation of {| 0 when Y = 0.1^> 
(filled circles) and for two cases where the gas and dust 
densities are not constant. The open circles and squares are 
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for non-uniform distributions of dust and gas. Tho open 
squares aro for a uniform dust to gas ratio with 
»(t) at n d (r) at o” T (or as functions of tho radius r, 
n ( r ) een d (r) at (1+ ft) , ft »• ^ r '* r u )<7 d r ) • Hie open circles 
aro for dust to gas ratio Increasing toward outer regions 
with n(T) «e e“ 2T and n d (T) * e" T (n(r) *(l+ft)' 2 and 
n d ( r ) « (1+ ft) ^) • 

Pig. 6, Variation of R, tho ratio of He to H volume emission measures; 
with Y 0 /V for various values of a (as labeled) and T^: 

^ « 0, no dust (dotted line), s 1,0 (solid linos) and 
T 1 " 2,0 (dashed lines). The other symbols aro as in fig. 

Hie open circles are connected to thoir corresponding lines 
for uniform nebulae by the vertical dashed linos. 

Fig. 7* Hie variation of R with the average value of ^ for various 
values of Y Q and T. , 
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Figure 1 
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Figure 2 
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Figure 4 
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Figure 5 
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Figure 6 
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